en 






Non- Hermit ian Hamiltonians in decoherence and equilibrium 

theory 

Mario Castagnino 
CONICET, lAFE (CONICET-UBA), 
IFIR and FCEN (UBA), Argentina. 

Sebastian Fortin 
CONICET, lAFE (CONICET-UBA) and FCEN (UBA), Argentina. 



O ; Abstract 

(N 

Vh I There are many formalisms to describe quantum decoherence. However, many of them give a 

^! 
■^ ' non general and ad hoc definition of "pointer basis" or "moving preferred basis", and this fact is 

a problem for the decoherence program. In this paper we will consider quantum systems under 

a general theoretical framework for decoherence and we will present a tentative definition of the 

*^' moving preferred basis. These ideas are implemented in a well-known open system model. The 

-)— > \ 

G . obtained decoherence and the relaxation times are defined and compared with those of the literature 



for the Lee- Friedrichs model. 
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I. INTRODUCTION 

From the appearance of quantum mechanics many attempts have been made to recover 
the laws of the classical mechanics through some classical limit. The more common scheme of 
this type includes quantum decoherence. This process eliminates the terms of interference of 
the density matrix, that are classically inadmissible, since they prevent the use of a classical 
(boolean) logic. In addition, decoherence gives a rule to select candidates for classical states. 

In this work the decoherence is considered an interaction process between an open quan- 
tum system and its environment. This process, called Environment-Induced Decoherence 
(EID) determines case by case which is the privileged basis, usually called moving preferred 
basis where decoherence takes place in a decoherence time to that is much smaller than the 
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relaxation time Ir and it defines certain observables that acquire classical characteristics. 
This is the orthodox position on the su 



case by case in several papers (see 



references 



fl 



ibject |l|. The moving preferred basis was introduced, 
I) in a non systematic way. On the other hand in 
3| and j^ Roland Omnes introduces a rigorous and almost general definition of 
the moving preferred basis based in a reasonable choice of relevant observables, and other 
physical considerations. Recently it has become evident that dissipation from system to en- 
vironment was not a necessary condition for decoherence |2| and the arrival to equilibrium of 
closed systems was also considered (|5[-[lJ]). Closed system will be discussed at large else- 
where. In this work we focus our attention on EID, which is a well-known theory, with well 
established experimental verifications, which makes unnecessary any further explanation. 

Non-unitary evolutions are essential to explain and study decoherence phenomena, quan- 
tum to classical limit, and final equilibrium. These phenomena appear in the evolution 
of quantum system, where decoherence time and relaxation time can be defined using non- 
unitary evolutions, poles theory, and non-Hermitian Hamiltonians. We will consider a closed 
system U and we will define two subsystems: 5, the "proper or open system" , and E, the 
environment. It is well-known that in this case the state of the proper system is obtained 
from the total density operator by tracing over the environmental degrees of freedom. If we 
consider the Hermitian Hamiltonian of a composed closed system U and the inner product 
of the evolved state with any observable we can make its analytical continuation, in the 
energy variable into the lower complex half-plane, and in general we will find poles. These 
poles are complex eigenvalues of the non Hermitian Hamiltonian Hefj that determines the 
system evolution. These complex eigenvalues define all the possible non-unitary decaying 
modes with characteristic decaying times proportional to the inverse of the imaginary part 
of the poles (see [l5|-[22]). From these characteristic times we can deduce the relaxation 
time is the largest characteristic time and it is related with the pole closest to the real axis. 
We can also deduce the decoherence time, that turns out to be a function of the imaginary 
part of the poles and the initial conditions of the system. Moreover, we will introduce a 
tentative definition of the moving preferred basis. All these definitions are considered in the 
Lee Friedrichs model. 



II. TOWARDS A DEFINITION OF THE MOVING PREFERRED BASIS. 

In this section we will try to introduce a very general theory for the moving preferred 
basis in the case of a general distribution of poles and for any relevant observable space 
Or. For this purpose it is necessary to consider the coordinates of observables and states 
in the Hamiltonian basis {\u))} (i.e. the functions 0{u),uj') and piuj.uj')) endowed with extra 
analytical properties in order to find the definition of a moving preferred basis. 

It is well-known that evolution towards equilibrium has two phases (there also is an initial 
non exponential Zeno-period which is irrelevant in this paper): 

i.- An exponential damping phase that can be described studying the analytical contin- 
uation of the Hamiltonian into the complex plane of the energy (see |l5l|-|22|). 

ii.- A final decaying inverse-polynomial in t~^ known as the long time evolution or Khalfin 



effect (see 



23| . |24j|). which is difficult to detect experimentally (see J25|). The power law 
decay for long times described by the Khalfin effect has no intrinsic parameter. It has no 
characteristic time scale. Khalfin period is the one where the decaying exponential modes 
are not more dominant and only inverse powers of time modes remain. We can consider that 
the time characteristic of this period is infinite (or very long). Instead of using the word 
"infinite", we will use Khalfin time scale. 

These two phases will play an important role in the definition of the moving preferred 
basis. They can be identified by the theory of analytical continuation of vectors, observables 
and states. To introduce the main equations we will make a short abstract of papers [l5| 



and |20|. 



A. Analytic continuations in the bra-ket language. 

We begin reviewing the analytical continuation for pure states. Let the Hamiltonian be 



H = Hq + V where the free Hamiltonian Hq satisfies ( see 15j or 20|) 



Hq\u) = uj\uj), {ojIHq = ijj{u\, < w < oo (1) 

and 

J= / duj\u){uj\, {u\u') = 6{u — u') (2) 

Jo 



Then 

poo 

Hq = I uj\uj) (uj\dbj (3) 

Jo 

and 



H = Hq + V = / u\uj){u\du+ du duj'V^^'\u){u'\ = / u\u^){u'^\du (4) 
Jo Jo Jo Jo 

where the \u}^) are the eigenvectors of H, that also satisfy eq. ([2]). The eigen vectors of H 
are given by the Lippmann-Schwinger equations (see |15J. eq. (12) and (13)) 

{^lj\u^) = {^\u) + {^\ ^^^^_^ V\uj), (a;+|v,) = (a;|^) + (a;|\/-— i— ^1^) (5) 

Let us now endow the function of u with adequate analytical properties (see |16|]). E.g. 
let us consider that the state \ip) (resp. {ip\) is such that it does not create poles in the 
complex extension of (wlv^) (resp. in {ip\uj)) and therefore this function is analytic in the 
whole complex plane. The physical meaning of this hypothesis is that if the system would be 
a non interacting one it would never reach equilibrium. Moreover we will consider that the 
complex extensions of function {u~^\ip) (resp. (iplu^)) is analytic but with just one simple 
pole at zq = uq — |7o, 70 > in the lower halfplane (resp. another pole Zq = Uo + |7o, 7o > 
on the upper halfplane) (see [lO| for details). Then in this paper, for the sake of simplicity 
we will always use a model with just one pole and an integral that corresponds to the Khalfin 
effect. Then we make an analytic continuation from the positive u axis to some curve F of 
the lower half-plane. 



Then (see (15|. eq. (29)) we can define 

(/olv?) = con0^2(,(u;'+|v9), (V^|/o) = {-27ii)cont^>^^^{u' - zo){tp\u^) 
ifz'l'f) = cont^>^^>{uj'^\ip), {^\fz') = cont^>^z{^\uj^), z' E T ,y \(p) {tp\ (6) 

and (see [l5|. eq. (31)) 

(^1/0) = cont^^^*{ilj\u+), {fo\ip) = {27ri)cont^>^^*{uj - Zo){u+\(p) 

(^1/2') = cont^^zi'ipl^^^), {fz\v) = cont^^z{uJ^\v), z eT,W \ip) (^| (7) 

where cont means analytic continuation. The tilde in (/o| is originated in the fact that in 
the complex extension there is no one-to-one correspondence between bra and kets [15| . 



Finally it can be proved that (see [l5| eq. 1.33 and J20| eq.82) 



H = zo\fo){fo\ + j z\fz){Mdz (8) 



That is a simple extension of the eigen-decomposition of H to the complex plane with the 
one-pole term and the integral term that produces the Khalfin effect. 

When it is possible to neglect the Khalfin term (i. e. for not extremely long times) the 
Hamiltonian reads where we have only a complex energy zq. 

H,ff = zo\fo){fo\ (9) 

This is the non Hermitian Hamiltonian that determines the evolution of the system far from 
the Khalfin time scale. 

B. Analytical continuation in the observables and states language. 

What we have said about the pure states and the Hamiltonian can be rephrased in the 



case of the states, observables, and the Liouvillian operator L (see a review in 27|). But 



we prefer to follow the line of [15| and keep the Hamiltonian framework and discuss the 
analytical continuation of (O)p(j), that we will also symbolize as (p(t)|0). In fact, we know 
that this scalar is the main character of the play so we will completely study its analytical 
properties. So let us call 

\uj) = \uj){ijj\, and la;,^;') = |u;)((^'| (10) 

Then a generic relevant observable is Or G Or 

Or = \Or)= f dLoO{cu) \co) + f dcu f dcu'0{uj, to') \cu, to') (11) 

and the generic states is 

pR = (prI = / dup{u){u\ + du du'p{u,u'){u,u'\ (12) 

where (see also (20| eq. (44) or [15|. eq. (45)). 

{(jj\(jj') = S{u — u') and {u,u'\ijj",ijj"') = 6{u — uj")5{uj' — uj"') (13) 

Then 

{^\0r) = 0{uj), {^'\0r) = 0{uj,uj') (14) 

We will consider the subject as general as possible, i.e. Or would be any observable such 
that Or G Or and pr any state pr G O'j^. In fact, in the next subsection we will only 
consider the generic mean value {pR{t)\OR) for two paradigmatic models below. Model 1 
with just one pole and the Khalfin effect and Model 2 with two poles. 



C. Model 1. One pole and the Khalfin term: 

We will use a formalism for states and observables which has been proposed by the 
Brussels school (led by Ilya Prigogine) in 26|. It can be proved (cf. ([l5|) eq. (67)) that the 



evolution equation of the mean value (p(t)|(9) is 

/'OO /"O 

(0«)p(i) = (p(t)|Ofi)= / p*{uj)0{uj)duj+ I p*itu,u')Oiu,u')e'^^-'''^'dudu' (15) 



oo poo 



'0 JO JO 

i.e. this mean value in the case V ^ reads 

(p(t)|0«)= /"rfc^(p(0)|$^)($:|Ofi)+ /"rfc^ /"dc^V(---')*(p(0)|<l>^^0(C'|Oij) (16) 

Where 0{oo) = {¥^\0n), 0{io,uj') = {^'\0r), p*{io) = (pr(0)|$,), p*{uj,io') = 
{Pr{0)\^ujuj'). These $ vectors are defined as 

1$^) = |a;+)(a;+|, I^.^O = |u;+)(u;+'|, (17) 

and 

(d^'|= [ d£[{uj+\6){6\uj'+) - 5(10 - 6)SiLo' - e)](i\ + I de ( de\uj+\e){e'\u'+)'{^'\ 

i%^\ = M (18) 

It should be emphasized that according to definitions flTUl) - flT^ , {io\ ^ {\uj^){uj^\) and 
{ujuj'\ 7^ (\uj^){ijj^'\) in contrast to the case of discrete spectra (see 15| for details). Then, if 



we endow the functions with analytical properties and there is just one pole zq in the lower 



halfplane, we can prove that (see [l5j eq. (70)) 

{p{t)\OR) = y"du;(p(0)|$^)($:|OH) + e*(^o-o)t(^(0)|$oo)(^|0«) 

+ /rfzV(^o-')*(p(o)|$o^,)($;^,|0^)+ f rfze^(-^°)*(p(0)|$o.)(5;;l|OR) 
Jv Jv 

+ f dz /dzV(^-^')*(p(0)|<l>,,0(5I'|OR) (19) 

Jr* JT 

where zq is the simple pole in the lower half-plane. |$^), {^z\i \^zz')i and {^zz'\ can be defined 
as in the case of eq. ([6j) and (jTj). The l^z), ("^zl, I'^zz'); ^^^1 ($22' | can be also defined as a 
simple generalization of the vectors |/o), (/o|, |/z), and (/z| ([1^. eq. (42)). 



Therefore we can conclude than the last four terms of equation (1191) vanish respectively 
with characteristic times 



12 2 

-;-; -;oo (20) 

7o 7o 7o 



Let us observe that: 



i. The vanishing of the second, third, and fourth therms of eq. (1191) is an exponential 
decaying corresponding to the first three terms of eq. f l20|) . This will also be the case 
in more complicated models with many poles. 

ii. The 00 in eq. (120|) means that the evolution of the last term of this equation corre- 
sponds to a polynomial in t~^ , i. e. to the Khalfin evolution. This is a very weak 



effect detected in 2006 25| . Therefore if there is a finite number of poles and the curve 



r, is below them, the contribution of the integral along F corresponds to the Khalfin 



effect. A closed system model for Khalfin effect can be found in 
an EID-like model in |[29j, section 5. 



section 6, and 



Then as we must have to ^ t^ and since from eq. fl20|l we have just two characteristic 
times 7(7^ and "00", the only possible choice is tn = 7^^"^ and t/j = 00. In fact, for times t 3> 
^D = lo^, eq. flT9|) reads 

{p{t)\OR) = j du{pm^u.){^u.\OR) + j dz j dz'e'^'-^'^\pm^zz'){^AOR) (21) 

since where t^ to = %^ the pole terms have vanished and we just have the Khalfin term. 
Let us now diagonalize p{t) of the last equation as 

pit) = jdzp,{t)m)m\ (22) 

where {|i(i))} is the moving eigenbasis oi p(t). Now let us define a state (pp(t)|, the preferred 
state, such that, for all times, it is 

{pp{t)\On)= /rfu;(p(0)|$J($:|O)+ / dz [ dz'e'^'-''^\pm<^,A{^AO) (23) 



r* 



So pp{t) is a state that evolves in a model with no poles and with just the Khalfin term. 
The functional (pp(t)| is defined by the inner product {pp{t)\Oji) as follows from the Riezs 
theorem^ . 

It is quite clear that 



^ All these formulas are confirnied by the coincidence of results with other methods: e.g. those used to 
study a '^'^^Pb{2d^/2) proton state in a Woods-Saxon potential (see [l5| Figure 3). 
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i. when t < tn, p{t) ^ pp(t) 

ii. when t — )■ to, p{t) — )■ pp{t) 

iii. when t ^ to, pif) = ppif) 

The eigen states of the pp{t) are those that we will choose for the moving decoherence 
basis. In fact, diagonalizing pp{t) we have 

pp{t) = Y.p,{t)m)m\ (24) 

i 

and when t ^ t£, = ^y^^ we have that p{t) — )■ pp{t) so from eqs. (l22l) and (12^ we see that 
the eigenbasis of p{t) and pp{t) also converge 

{W))]^{\m)} (25) 

Namely the basis {|^(t))} converges to {\j{t))} and therefore p{t) becomes diagonal in 
{|j(t))}. Thus {|j(t))} is our definition for the moving preferred basis for this case. Since 
pit) becomes diagonal in the just defined preferred basis {|j(t))} when t —> Id and t^ = 7^^"^ 
is really the definition of the decoherence time. In this model the relaxation time tf> corre- 
sponds with the Khalfin term, i.e. an extremely long time, so that 

tn < tR (26) 

D. Model 2: Two poles without the Khalfin term. 



The Khalfin term is so small (see [25||) that it can be neglected in most of the experimental 
cases. Then we can eliminate the Khalfin term since it corresponds to extremely long time. 
In this case the Hamiltonian becomes non-Hermitian as in eqs. (jH]) and (IM|) . So let us 
consider the case of two poles zq and zi (and no relevant Khalfin term) where eq. flT^ reads: 



ipit)\On) = j rfu;(p(0)|$^)($JOK) + e*(^o-^")*(p(0)|$oo)($oo|0^) + 

+ e^(^i*-^«)*(p(0)|<l>io)($^|OR) + e*(^o-i)*(p(o)|$oi)($;;;|0^) 

+ e*(%*-o)t(p(0)|<l>n)($;;|0«) (27) 



where zq = Uq — ^■Jq, ■Jq > , Zi = Ui — ^'ji, 71 > 0, and we will also consider that 70 ^ 71 



(see 30| section 3, for details). Then the characteristic times (I2CT|) now read 

1111 

— ; = ^ — (28 

7o 7i + 7o 7i + 7o 7i 

Then we must choose ta = 7f ^ and t^ = 'Jq'^. Now for times t ^ tn = 'Ji'^, eq. f l2T|) reads 



ipit)\OR) = j da;(p(0)|<l>J($JO^) + e^(^o*-o)*(^(o)|$oo)($oo|OR) (29) 

and we can define a state (pp(t)| such that, it would be 

ippit)\On) = I duipiO)\<^^){^^\On) + e'^^o-^°^\piO)\%o)i%o\OR) (30) 



for all times. Repeating the reasoning from eqs. fl2T|) to fl25|l we can see that, diagonalizing 
this last equation, as in eq. (12^ . we obtain the moving preferred basis. Then in this case 
we see that the relaxation is obtained by an exponential damping (not a Khalfin term) and 

ifl = — > tD = — (31) 

7o 7i 

Again, in this case when t ^ tr, = ^^^ we have that pnit) — )■ pp(t), and we can conclude 
that the eigenbasis of p{t) and pp{t) also converge as in eq. fl2^ . Namely p{t) becomes 
diagonal in the moving preferred basis in a time to- 

E. The general case 

Let us now consider the general case of a system with A^ + 1 poles at Zi = u'^ — i'ji. In 
this case it is easy to see that eq. (ITTI) (with no Khalfin term) becomes: 

N 

{p{t)\Op) = {p,\Or) + Y, <t) exp (-7it) = ipR.\Op) + fit) (32) 

j=0 

where (p*|Or) is the final equilibrium value of (p(t)|Oi?) and the ai{t) are oscillating 
functions. In the most general case the Zi will be placed either at random or following 
some laws. Anyhow in both cases they can be ordered as^ 

70 < 71 < 72 < ■•• (33) 



^ For simplicity we will only consider the case 70 ^ 71 ^ 72 ^ ...Other special cases will be considered 
elsewhere. 
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Then if 70 <^ 71 it is quite clear that the relaxation time is 

tR=- (34) 

7o 

So the relaxation time is defined with no ambiguity. Let us now consider the decoherence 
time. Really each pole zi defines a decaying mode with characteristic time ti = 7"^. More- 
over the poles contain the essence of the decaying phenomenon and the definition of the 
decoherence time depends on their distribution and other data like the initial conditions. 
In fact, the initial conditions seam essential for the definition of to- To introduce these 
conditions, let us define: 

TV N 

fit) = Y, a^{t)e-''^\ fit) = Y, a'Me-^^' - a.(t)7.e-^^* (35) 

i=0 i=0 

SO at t = we can write the initial conditions as 

N N N 

/(O) = Y «^(0)' /'(O) = E ^^(0) - E «^(0)^^ (36) 

i=0 1=0 i=0 

Let us call f{t) = const, exp g{t) ~ expg{t), and let us make a Taylor expansion of g{t) as 

g{t)=g{0) + g\0)t+y'{0)t' + ... (37) 

So let us postulate the reasonable hypothesis that the decoherence time is to '^tf>. Then, 
in the period before decoherence that we are considering , precisely t < to "^ tn, we have 
^ ^ 1. With this condition we have the approximation: 

g{t)=g{0) + g'iO)t (38) 

where 

^(O) = log/(0), ^0) = ^ (39) 

These equations contain the initial conditions. Then in this approximation: 

/(t) = e^(o)eV(o) ^ j(o) exp f ^f^t | exp ( -^^^^fllt] (40) 



So we define 



.2^1=0 '^« / \ l^i=o'^i 



a.(t) = /(0)exp|^f^t and 7e// = ^^f^ (41) 

11 



And (HUl) becomes 

/(t) = a,(t)exp(-7e//t) (42) 

The decoherence time is 

tD = ^ (43) 

7e// 

Then 7e// and to are both functions of the initial conditions. We will see that this to 
coincides with the one of the Omnes example in the next subsection. 

Let us now consider the definition of the moving preferred basis. It is clear that, for the 
time t 3> t_D5 the modes with characteristic times tj < tu (i.e. ■ji > leff), that we will call 
the fast modes, have become negligible in eq. (132|) . Then we can define the functional (pp(t)| 
as 

M 

{pp{t)\OR) = (p,|Or) + 5^a,(t)exp(-7,t) (44) 

i=0 

where the sum in this equation only contains the M < N poles such that 7, < 7e//, where 
the 7j correspond to the slow modes. This is our adiahatic choice since we have selected 
the slow modes of decaying to define pp{t) and rejected the fast modes. Our adiabatic 
choice corresponds to keep the slow modes and disregard the fast ones. Thus, for us the 
robust modes are the slow modes since they are "the less affected by the interaction with 
the environment" , that creates the poles, if compared with the fast modes, and it is usual 
to say that these robust modes are those that define the moving preferred basis. In fact: 

i.- If the Hamiltonian would only be Hq (cf. eq. (|3])) there would not be poles (and this 
is the usual case in the literature). But the complex extension of the complete Hamiltonian 
H (cf. eq. (3)) certainly has poles. Therefore the poles are created by the interaction 
Hamiltonian V . 

ii.- Thus the slow modes and the fast ones are defined by these poles, and in the case we 
are considering, i.e. EID, the poles are defined by the interaction with the environment. 

iii.- Then it is reasonable to call robust the slow modes, since the environment interaction 
has smaller influence in these poles, and we conclude that these are the modes that define 
the moving preferred basis. 

This is our definition of robustness. Analogously, if we compute the linear entropy we 
will have a slower variation of this entropy, if we only consider the slow modes, than if we 
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consider all the modes (including the fast ones). This would be our minimization of the 
linear entropy: the moving preferred basis evolution only contains the slow modes. 

Moreover, when t ^ tr, the motions produced by the fast modes, such that 7^ > 7e//, 
namely those with motions faster than the one of the evolution of eq. fHTl) . are no more 
relevant for p(t), and pp{t) — > p{t). Then we diagonalize pp{t) and we obtain the moving 
preferred basis {\i{t))}. The only influence in the evolution of pp{t) is given the poles such 
that 7j < 7e//- When t — > t^, {K(^))} ~^ {b(^))} the eigenbasis of p{t) where < t < 00. 
This {\j{t))} is our candidate for a general definition of moving preferred basis. 

III. THE OMNES OR LEE-FRIEDRICHS MODEL. 

Our more complete and sirnplest example of decoherence in open systems is the Omnes 
"pendulum" (i. e. oscillator jSli) in a bath of oscillators, that we will compare with the 
poles theory in the following subsections. In fact the Omnes model could be considered a 
poles model if we retain the poles and neglect the Khalfin term. Moreover in the Omnes 
philosophy the moving preferred basis must be related to some "collective variables" in such 
a way that they would be experimentally accessible. In this case this variable is the center of 
mass of the pendulum, i. e. the mean value of the position of a coherent state. In [31l| page 
285 a one dimensional "pendulum" (the system) in a bath of oscillators (the environment) 
is considered. Then the Hamiltonian reads 

H = ua^a + Y^ ujkblbk + ^(Afca"^6fc + Xlabl) (45) 

k k 

where a'^{a) is the creation (annihilation) operator for the system, h\{hk) are the creation 
(annihilation) operators for each mode of the environment, u and Uk are the energies of the 
system and of each mode of the environment and A^ are the interaction coefficients. 
Then let us consider a state 

\^{t)) = a\a,{t)) n |/3fci(t)) + 6||«2(t)) n \h2{t)) (46) 

k k 

where |ai(0)), |«2(0)) are coherent states for the "system" corresponding to the operator a\ 
with center in Xi(0) and ^2(0) respectively, and |/3a;i(0)), |/3a:2(0)) are a coherent state for 
the environment corresponding to the operator h\. Let the initial conditions be 

IV'(O)) = a|ai(0) {hm = 0}) + 6|a2(t), {/3fc2(0) = 0) (47) 
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Moreover Onines shows that, under reasonable hypotheses and approximations the relaxation 
time of the system is 

tR = 1/7 (48) 

where 

7 = 7r /" n{v')dv'\l,5{co - v') (49) 



where n{v')dv' = dk. On the other hand, the decoherence time of the system is (see J3l|, 
pp. 289-291) 

in - -^tR (50) 

where Lq =|xi(0) — X2(0)|. In the next subsection, we will attempt to recover these results 
using the polar technique. 

A. The characteristic times from the polar technique. 



A particular important model can be studied, like the one in [211], with the Hamiltonian 
H = cuQa''a+ / cokbl^bjid'k + / Afc(a^6k + 6|^a)(ik (51) 



i.e. a continuous version of (l45l) . In this continuous version we are forced to endow the 
scalar (p(t)|Oij) with some analyticity conditions. Precisely function A^ (where k = uJk = |k|) 
is chosen in such a way that 

/ N f dkXl 
'n±[oJk) =uJk-uJo- / —-7: (52) 

does not vanish when k G M^, and its analytic extension ^+(-2) in the lower half plane only 
has a simple pole at zq. This fact will have influence on the poles of (p(t)|0/j) as in the 
last section and we know that the study of {p{t)\OR) is the essential way to understand the 
whole problem. 

The Hamiltonian flSTl) is sometimes called the Lee-Friedrichs Hamiltonian and it is char- 
acterized by the fact that it contains different number of modes sector (number of particle 
sectors in QFT). In fact, a) and h\^ are creation operators that allow to define these numbers 
of mode sectors, e. g. the one mode sector will contain states like a^\Q) and 6^10) (where 
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o|0) = &k|0) = 0). Then the action of exp {—Ht) (or simple the one of H) will conserve the 
number of modes of this sector in just one mode, since in flHT]) all the annihilation operators 
are preceded by a creation operator. This is also the case for the n— mode sector. 

1. The Friedrichs model and the relaxation time 

In the case of the one mode sector this model is the so called Friedrichs- Fano- Anderson 



or Friedrichs model. For a complete discussion on this model see 32|. The Hamiltonian of 
the Friedrichs model is 

Hf = uo |1) (1| + Uk \uj) {uj\ duj+ I {\{w) \uj) (1| + \*{u) |1) (wl) duj (53) 

(this Hamiltonian which is similar to the one of eq. dl]), is expressed just in the variable u, 
the one that will be analytically continued). As a consequence of the analyticity condition, 
mentioned above, this simple Friedrichs model only shows one resonance. In fact, this 
resonance is produced in zq. In paper t32i] we can see that the poles we compute here are 
the same as the poles of the Green's function. Let Hp be the Hamiltonian of the complex 
extended Friedrichs model, i.e. the Hamiltonian of eq. (|8]), then^ : 

Hf\zo) = zo\zo), Hf\z) = z\z) (54) 

where zq = ujq + 5uJo — ijo = u'q — i'jo is the only pole and z G F corresponds to the integral 
term and to the Khalfin effect. 

The Lee-Friedrichs model, describing the interaction between a quantum oscillator and 
a scalar field, is extensively analyzed in the literature. Generally, this model is studied by 
analyzing the one excited mode sector, i.e. the Friedrichs model. Then, if we compute the 
pole, of this last model, up to the second order in Xj. we obtain that 

zo = uo+ ^ 55 

So the pole (that will corresponds to the one closest to the real axis in the Lee-Friedrichs 



model) can be computed (see |33| eq. (42)). These results coincide (mutatis mutandis) with 
the one of Omnes book ^] page 288, for the relaxation time. In fact: 

^ p( ^ ] - in5{LOo - tu') (56) 



Uq — u' + iO \Uo — u' 



Only symbolically, since the poles really belong to the scalar (p(i)|0), as in the last section. 
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where P symbolizes the "principal part", so 

dk'Xl 



zo = uJo + P / in / dk'X^,6{uo — ujk) (57) 

J UJq — UJk J 

Then if dh = n{u)duj we have 

S^^p f n{u')du'\l ^ ^ ^ [n{uj')dio'Xl,6{uJo - uj') (58) 

J LOq — LO J 

where 5uj is a shift and 7 a damping coefficient, then the system would arrive to a state of 
equilibrium, namely the results of [31] page 288, and the one contained in eq. ( H9|) yields: 

Zq = {uq + 5u) - 27 = Wq - i7 (59) 

So the Omnes result for the relaxation time coincides, as we have already said, with the one 
obtained by the pole theory, precisely (see (148|) ) 



tR = - (60) 

7 

2. Other poles of the Lee-Friedrichs model. 

Let us now consider the Lee-Friedrichs Hamiltonian f lSTj) for the many modes sectors. 
Then, as an example for the three mode sector (with just the unique pole zq and Zi, Z2, or 
Z3 "real continuous eigenvalues" transported to the curve F) we have: 

H\Za, Zb, Zc) = {Za + ^6 + Zc)\Za, Zb, Zc) (61) 

where {za + Zb + z^) is the eigenvalue. Then 2:1,^2; ^3 € F is the Khalfin terms (i.e. they 
belong to the complex contour on the lower complex energy plane), and let Zo be the pole 
of one particle sector. So in the real complex plane the spectrum of H contains 

1.- Eigenvalues {zi + Z2 + -23) with three points of the curve F. 

2.- Eigenvalues (zi + Z2 + zq), {zi + zq + z^) and {zq + Z2 + z^), with two points of the 
curve F and the pole Zq. 

3.- Eigenvalues {zi + ^o + ^o); (-^o + -^2 + -^o) and (zq + ^o + ^3), with a pole at 2zo, and 
one point of the curve F. 

4.- Eigenvalue (2:0 + zq + zq), with a pole at 3zo. 

These values appear in expression of the mean value as ~ e~*"fi"* (like in eq. (19) second 
term of the l.h.s.) or as ~ J-^e~^'i^^f{z)dz (like in eq. (19) three last terms of the l.h.s.). 
Then we have that the four cases above become: 
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!•- It It It ^ ' ' ^ ^ */(^i, ^2, Z3)dzidz2dz3 

2.- J-^J^e'"^ h ^f[zi,Z2)dzidz2, and the same for the combinations {zi + zq + Z3), 

{Zq + Z2 + Z^) 

3.- Jp e~* ft ^f{zi)dzi, and the same for the combinations (^2 — 22;o), (,23 + 2zq) 



, 3zn 



4.- e-^* 

Then if we neglect the Khalfin we just have the point 4. 

Of course in the general case 3 — > n we would have e~* « (for the point n) , plus many 
integrals on the curve F (for the points 1,2, ...n — 1) corresponding to Khalfin terms Then 
if we neglect the integrals that produce the Khalfin effect, since this effect corresponds to 
extremely long times, the F term disappears and we simply have a pole at Zn = uzq. This 
elimination (in the case of just one pole zq) introduces in the model a structure of a complex 
oscillator. Then we can introduce a non Hermitian effective Hamiltonian 



Heff = zo ( ajao + - j = zq (No + -j 



(62) 



where Oq and a^ are creation and annihilation operators and A^o = % ao is the number of 
modes operator and 

NqIu) = n\n) (63) 

In the case of large n, H^ff becomes extremely close to 

Heff = zoalao = zqNq (64) 



Moreover we can call 



nzQ = n{ujo - i'jo) (65) 



and we will find the evolutions 

exp{—iHefft)\n) = exp{—inzQt)\n) = exp{—iznt)\n) (66) 

So, in this approximation, the effective Lee-Friedrichs Hamiltonian H^ff simply is a (non 
Hermitian) version of H with just damping terms. We will below use this structure. 
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3. The initial conditions 

As an initial conditions, |ai(0)), |a2(0)), it is possible to choose any linear combination 
of the elements {\n)} where n = 0, 1..., oo. So we can choose coherent states 

l«.(o)) = e-^ Y. -^rr- 1^^ (^^) 



n=0 

Then let us choose the initial conditions as the sum of two coherent states, namely: 

|$(0)) = a|ai(0)) + 6|a2(0)) (68) 

Thus the initial state operator is: 
p(0) = |a|' |ai(0)) (ai(0)| + |6|' 1^2(0)) (a2(0)|+a6* |ai(0)) (a2(0)| +a*6 1^2(0)) (ai(0)| (69) 



We choose the two Gaussian ( 1^ with center at ]5i,2(0) = 0, (see [3l| eq. (7.15) page 284) 
and 

TTIOJ TflOJ 

ai(0) = ^^xi(0) and ^2(0) = ^^0:2(0) (70) 

y2muj y2muo 

So ai(0) and 0:2(0) are real numbers. With a change of coordinates we can choose xi(0) 
and ^2(0) without loss of generality. So we can consider that the ai(0) and 0:2(0) are 
both positive. For this reason we will interchange Oj(0) and |oj(0)| below. With no lost of 
generality we can choose 

TflUJ 

Oi(0) = and 02(0) = Lo (71) 

\/2muj 

The macroscopic case 

It is easy to prove that for macroscopic initial conditions, i.e. when the peaks of the 
two Gaussians are far from each other, that is to say |oi(0) — 02(0) | — )■ 00, the states 
{|oi(0)) , |o2(0))} are quasi-orthogonal basis 

(ai(0)|o2(0)) - (O2(0)|oi(0)) - 6- '"^'°^"^°^'°"' - (72) 

Then, the macroscopic condition for the initial conditions is |oi(0) — O2(0)| ^ 1. So we 
have |oi(0) - 02(0)] = 02(0) > 1, 

"" i„ » 1 (73) 



\j2mijj 
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4- Components of the non- diagonal part of the state 

Therefore the evolved state is 

p{t) = \af |ai(t)) (ai(t)l + \bf ^t)) (asWI + ab* \ai{t)) {a2{t)\ + a*b Mt)) {ai{t)\ (74) 

Let us not consider the non-diagonal part of p(t), p^^^\t) in the basis of the initial conditions 
{|a;i(0)) , |a;2(0))}. Then we have 

p(^^)(t) = pSr^(t) MO)) {a2m+pT\t) l«2(0)) («i(0)| (75) 

Since the basis {|q:i(0)) , |a2(0))} is quasi-orthogonal, from eq. fITil) we have 

P^''\t) = ab* (ai(0)|ai(t)) (a2(t)|a,(0)) + a*b {ai{0)\a2{t)) (ai(t)|a,(0)) (76) 

If we consider the evolution given by the non Hermitian Hamiltonian 

|«,(t)) = e-^^^//* |a.(0)) = e-^^ f] i^l^e-^"' \n) (77) 



n=0 * 

we can compute these products {ai{0)\aj(t)) and we can replace them in fl76l) to obtain 

pjr^(t)-a6*e^"^(°)l'H'"'0 
p(7^)(t)-a*6e-l"^(°)l'(i-^"'^"') (78) 

5. Decoherence time 

Since the contributions 7„ of individual poles Zn do not appear explicitly in the equation 
(ITS]) , we may think that such poles are not involved in the outcome. 

However, if we express the exponential of fITS]) as its Taylor series, we have 

pSr^(i) = a6*e-l-(°)l^ y M£lQ!e<* = ab* V c„(t)e-^"* 

n=0 * n=0 

p(f^)(t) - a*6e-l"^(o)l' f; M^e-"* = a*6f;<(t)e-^"* (79) 



where 



n=0 n=0 



|a2(0)P (10^2(0)!')" ^,^„, 



c„(t) = e-l"2^"^l ^' ^ ^^' ' e'""* (80) 
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and these equations show that all the diagonal terms vanish when t — )■ oo showing that there 
is decoherence. Now we would like to know the decoherence time, then we must find 7e//- 
So we analyze the decay of 

pr'wf=psr'w(psrw)*=i«i.i^»"' (sd 

where 

9it) = In ( f £ c^{t)e'^A f £ c:{t)e-A j (82) 

Let us now expand e^^*-* as: 

(83) 



,9W _ ^ff(0)+5'(0)t+|g"(0)t2 + . 



As the decoherence time is a very short one to ^ t^ let us neglect it from the quadratic 
term. Now, from eq. fl82|) we have that 

9(0) = 

/(0) = -2e-l«"»pf: 1^510!^,. = ,,„ (84) 

71=0 

then from (I5T1) and (IMl) . we have 

p!J^)(t)|-|a6|exp(-7e//t) (85) 

This is precisely the interpolation that corresponds to eq. (1421) . Now we have the decoherence 

time 

1 2 1 

In fact this «B' turns out to be a weighted average of the .maginary part of the poles .„. 
The same time was found by Omnes in [SlJ (or in eq. ([SD]) of this paper) and corresponds 
to the definition 0431) of the last section. So in fact, we have found the same result. Also 



m 



3l| the result for to is only valid for small t as in the last section. So the coincidence of 
both formalisms is proved. 

Let us now consider the mathematical definition of moving preferred basis. The basis 
{\ai{t)) , \a2{t))} = {\ai(t))}, is orthonormal when L — > oo and the reasoning below is done 
under this condition. Then let us diagonalize pp(0) (always when L — )■ oo) as 



pp{t) = $^p.(t)K(t))(z(t)| =pi(t)|l(t))(l(t)|+p2(t)|2(t))(2(t)| (87) 



i=l,2 
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where {\i{t))} is our orthogonal moving pointer basis. But for times t':^t£, pR{t) = p{tY^^ 
we have 

pp{t) = p^''\t) = |a|2 |«i(t)) {am + \h\^ |«2(t)) (a2(t)| (88) 

Then, since a hnear orthonormal decomposition is unique we find the moving pointer basis 

\l(t)) = \a,{t)). |2(tj) = l«2(t)) (89) 

Then Omnes basis coincides with {|l(t)), |2(t))}, but this also is our moving preferred basis 
since it evolves under the slow motion pole evolution. So Omnes basis and our basis coincide 
(always when L — t- oo). 

So we have proved that all the characters of the Omnes model: t/j, Id, and the moving 
preferred basis, coincides with our definitions of the last section. 

IV. CONCLUSIONS 

In this work we have: 

i.- Discussed a general scheme for decoherence, that in principle could be used in many 
examples. 

ii.- We have given a quite general definition of moving preferred basis {|j(t))}, and of 
relaxation and decoherence times for a generic system. 

iii.- We have proved that our definitions coincide with those of the Omnes model. 

We hope that these general results will produce some light in the general problem of 
decoherence. 

The Omnes formalism, of references |3|, [^, and 3l| contain the most general definition 



of moving preferred basis of the literature on the subject. Our basis has another conceptual 
frame: the catalogue of decaying modes in the non-unitary evolution of a quantum system. 
But since the Omnes formalism is the best available it would be very important for us to 
show, in the future, the coincidence of both formalisms, as we have at least done for one 
model in this paper. 

Of course we are fully aware that, to prove our proposal, more examples must be added, 
as we will do elsewhere. But we also believe that we have a good point of depart. In fact, 
probably the coincidences that we have found in the Omnes model could be a general feature 
of the decoherence phenomenon. Essentially because, being the poles catalogue the one that 
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contains all the possible decaying modes of the non unitary evolutions, since relaxation and 
decoherence are non-unitary evolutions, necessarily they must be contained in this catalogue. 
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